Quantum-limited position measurements of a dark matter-wave soliton 
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We show that the position of a dark matter-wave soliton can be determined with a precision that 
scales with the atomic density as n -3 ^ 4 . This surpasses the standard shot-noise detection limit for 
independent particles, without use of squeezing and entanglement, and it suggests that interactions 
among particles may present new advantages in high-precision metrology. We also take into account 
quantum density fluctuations due to phonon and Goldstone modes and we show that they, somewhat 
unexpectedly, actually improve the resolution. This happens because the fluctuations depend on 
the soliton position and make a larger amount of information available. 
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I. INTRODUCTION 

Degenerate quantum gases of atoms and molecules can 
be trapped in space and their interactions can be con- 
trolled by laser and magnetic fields. The spatial quan- 
tum state in these gases can be probed in ways un- 
paralleled in other many-body systems. A number of 
experiments have thus been devoted to exploring their 
coherence properties by the observation of interference 
patterns [lj and topological excitations such as solitons 
and vortices [1, H, [j] • For application of ultra-cold gases 
in high-precision metrology, it is a central question how 
to extract the maximum information from the accessi- 
ble measurements. We address the case of spatial den- 
sity measurements on a Bose-Einstein condensate (BEC) 
with a dark soliton dip, which is a minimum in the atomic 
density profile. Dark solitons have been produced by 
slowing down a light pulse using electromagnetically in- 
duced transparency (EIT) [J] , and the precision probing 
of the position of such solitons may be used to measure 
detailed properties of EIT. 

It has also been proposed, with a similar setting as 
in Ref. to use a time-dependent atomic trap guide 
as a matter-wave interferometer where, after splitting 
and recombining of a BEC, an oscillating dark soliton 
is formed @ . A high-precision determination of the soli- 
ton position is crucial here since the amplitude of the 
soliton oscillation depends on the interferometric phase. 
The signal-to-noise ratio is in many interferometric ex- 
periments given by the standard shot noise, leading to 
a phase resolution that scales with X/yN, where N is 
the number of detection events. The normal interfer- 
ence fringes in a matter-wave interferometer are located 
in accordance with the interferometer phase and yield the 
same dependence on the number density. Special quan- 
tum features such as entanglement and squeezing imply 
that this resolution limit is not a fundamental one Q ■ In 



principle, the Heisenberg uncertainty relation provides 
the optimum sensitivity with a scaling as l/N [8j, and 
for particular entangled states the interferometric phase 
error scales as l/iV 3 / 4 Q. Our work will show that the 
position of a dark soliton can be determined with a pre- 
cision scaling with the same power law. This result, how- 
ever, does not rely on squeezing or entanglement. It can 
be qualitatively understood from counting shot noise, the 
slope of the mean density versus position, and the 1/ ' ^fn 
dependence of the soliton width on atomic density due to 
interactions. Our analysis determines the Fisher informa- 
tion available in the spatial recording of the atomic den- 
sity and the resulting Cramer-Rao lower bound (CRB) 
[1 01 ] on the determination of the position of the soliton 
dip. 



II. SPATIAL IMAGE OF A DARK SOLITON 

The experimental situation we have in mind is sketched 
in Fig. [T]Ja). A quasi one-dimensional (ID) BEC is im- 
aged by counting the number of atoms present in sub- 
intervals, referred to as pixels in the following, spanning 
the extent of the system. For simplicity we imagine the 
probing to occur with unit efficiency, and we note that 
the recorded data will fluctuate around the quantum- 
mechanical expectation value of the number of atoms 
present in each pixel. 

A dark soliton is a stationary excited state solution of 
the Gross-Pitaevskii equation (GPE) [Tlj . 
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where $(cc) is the order parameter, m is the atomic mass, 
g is the effective ID interaction strength, and fi is the 
chemical potential. For finite systems <& is normalized to 
the number of particles in the gas, and is proportional 
to the single-particle wave function, self-consistently oc- 
cupied by all the atoms in a Hartree-Fock product state. 
When g > (repulsive interactions) the dark soliton in 
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Figure 1: (Color online), (a) Dark soliton density distribu- 
tion (shaded region) and typical measurement outcome with 
number fluctuations (red thick curve), (b) Rescaled Fisher in- 
formation F£ 2 versus pixel size for a mean-field solution with 
Poisson counting statistics (solid curve) and with a Gaus- 
sian continuous approximation (dashed curve) . The optimum 
sensitivity for locating the soliton position is Ag = F^ 1 ^ 2 ; 
£ oc n _1//2 is the healing length, (c) Rescaled Fisher informa- 
tion versus atomic density: Poissonian distribution (circles), 
Bogoliubov and Gaussian approximations (diamonds). Con- 
tinuous thick (black) line represents the behaviour of Eq. Q, 
whereas the dashed (green) line shows the information ex- 
tracted from the signal-to-noise ratio for a simple gain func- 
tion g(x,n) (see text). 



the homogeneous case is given by [TTjj 

q) ~ \/«tanh [k(x — q)] , (2) 

where n = l/(\/2£;) is the reciprocal width of the soliton 
profile, n is the linear density of the condensate, and 
£ = h/y/2mgn is the healing length. The position of 
the soliton is described by the parameter q. 

When dealing with not truly ID Bose gases, even in 
the dimensionality cross-over [12|], care has to be taken. 



Exact ID solitons are stable because of the integrability 
of the ID GPE, while solitons in 3D elongated conden- 
sates are thermodynamically unstable because of scat- 
tering with thermal excitations [HI]. We consider our ID 
model to be a good description for the soliton dip un- 
certainty, especially when a measurement is performed 
shortly after its preparation [l4| ■ We note also that part 
of our analysis can be readily generalized to the full 3D 
description and also to effective ID descriptions with a 
suitably modified GPE [H, 

The number n s of atoms in pixel s will on average be 
equal to the integral of the density p(x;q) = \&(x;q)\ 2 
over the pixel area Ax: n s — f* B+Ax dx p(x; q). Due 
to fluctuations around the quantum-mechanical expecta- 
tion value, the counting data will be noisy, and we need 
a theory that provides the best estimator for q. The 
Cramer-Rao bound provides a lower limit on the min- 
imum uncertainty of a parameter X, determined from 
measured values of a quantity 77, whose conditional prob- 
ability distribution is given by p(r)\X). The variance of 
an unbiased estimator of X is limited from below by 
min[Var(A)] = l/F(X), where 



F(X) = - [drip^X 



d 2 lnp(rj\X) 
dX 2 



(3) 



is the Fisher information [lfj. Let us assume that q is 
known (by a coarse scale estimate) to within an inter- 
val for which a first-order series expansion yj p(x; q) ~ 
yj p(x; 0) + qf(x) is precise enough. A mean- field con- 
densate, whether in a total number eigenstate or a co- 
herent state, displays Poissonian counting statistics in 
each pixel, when registered in a large number of them. 
An estimator that reaches the CRB is found by weight- 
ing the measured data n s with a weighting function 
g(x) oc f(x)/p(x;0) [l7[. We thus have a procedure by 
which we can estimate the soliton position optimally. 

To identify the CRB we need to evaluate the Fisher 
information ([3|), where rj is the entire set of (discrete) 
variables {n s } with a probability distribution that de- 
pends on the soliton position q. The Fisher information 
can be determined analytically, as done in Ref. [171 ] for a 
coherent state of light populating a mode with a given 
spatial density. For our dark soliton @ we get 



F = 4 dx 



dy/p(x;q) 



3ft 



(4) 



This dependence on the density is shown in Fig. [ljc) 
(continuous thick line) . The main result is that the small- 
est value of q that can be distinguished scales as n -3 / 4 
with the atomic density. This is a faster convergence 
with n than the usual shot- noise scaling (~ n -1 / 2 ). We 
emphasize that this enhancement does not require any 
squeezed or otherwise entangled multiatom state; it fol- 
lows from a simple (coherent state) macroscopic wave 
function. Note also that we obtained Eq. (0J by assuming 
no correlations between different pixels and taking the 



3 



limit of infinitely small pixel areas. Fig. [ljb) (thick solid 
line) shows the dependence of the rescaled Fisher infor- 
mation F £ 2 on the pixel size (placing a soliton dip at the 
border between two pixels). 

The dark soliton owes its existence to interactions, but 
interactions also cause deviations of the exact quantum 
many-body state from the simple mean-field solution. 
These deviations cause quantum depletion of the conden- 
sate and fluctuations of the soliton dip position, which 
cannot simply be set to a classical value. We analyze 
these effects by Bogoliubov theory, and show that, con- 
trary to what may be expected, quantum fluctuations 
actually increase the optimum sensitivity of the soliton 
position measurement. 

First, we note that the case of Poissonian count- 
ing statistics is tractable because the integral (sum) 
in Eq.((3]) over all possible measurement outcomes can 
be performed analytically. Poisson statistics is fully 
parametrized by the mean values of the observables, and 
hence it does not offer the possibility of investigating the 
role of extra noise or correlations in the gas. We therefore 
establish a formalism in which it is possible to deviate 
from Poisson statistics and in which a good approxima- 
tion to the Fisher information ([3]) can be computed effi- 
ciently As a natural choice we consider a multi-variate 
Gaussian approximation for the distribution function of 
the atomic density, namely 



(2tt)- a W 2 

p{nk) = -vmr exp 



— (n-n^C-^n-n) 



.(5) 



Here the column vectors n and n collect the detected 
stochastic atom number variables and their expectation 
values, and M px is the total number of pixels. The corre- 
lations between atom numbers on the pixels are described 
by the covariance matrix C. 

Inserting ([5]) in Q gives the Fisher information as 



F = - 
2 



1 j d 2 det{C)/dq 2 



det(C) 



ddet{C)/dg 
det(C) 



(6) 



E 



When Eq.© is evaluated for a continuous Gaussian dis 
tribution with the same covariance as a discrete Pois- 
son distribution, C s j = Var(n s ) 5, 
Fisher information 



SJ = n s S S j, we get the 



1 / dn t 



1 ^ 1 / 9n s 

9 i2 



(7) 



The first sum is equivalent to the Poisson result, while 
the second sum is a consequence of our Gaussian approx- 
imation. In the limit n s 3> 1, where the distributions 
are similar, this latter term becomes a negligible correc- 
tion. For infinitely small pixels, however, the Gaussian 
and Poissonian distributions differ and the second term 



in Eq.([7|) diverges as shown in Fig. [ljb). This diver- 
gence is linked to the nonphysical, noninteger, and neg- 
ative values of n allowed by the Gaussian distribution. 
In the analysis of real experiments, however, this is not 
relevant, since detectors will typically not be chosen so 
small that only one or zero atoms are recorded in some 
pixels. 



III. BOGOLIUBOV-DE GENNES THEORY 

In order to determine the quantum fluctuations of the 
atomic density, we have diagonalized the Bogoliubov-de 
Gennes operator 



C 



BdG 



(H G ?+g\<$>\ 2 
I -g<S>* 2 



g$ 2 
-H GP -g\<S>\ 2 



(8) 



Here, the wave function $ = <I>(x; q) is the stationary 
mean- field solution given by Eq.Q. Following the ap- 
proach in Rcfs.fTil. Il8j |. we find that the phonon mode 
functions (uk,Vk) = (W^ , Wj~) of £BdG m a box with 
periodic boundary conditions can be written as 



K 



/?k (fc/2 + i k tanh[/t (x - 



?)] + 
«)])}, 



(9) 



where is a normalisation constant [l9|, (3^ = (k/ k) 2 ± 
2e k /(gn), and e k = Hc\k\ y/l + k 2 /(An 2 ). 

Before discussing the contribution of the phonon 
modes, we note that Cbag m addition has two gapless 
modes. They arc Goldstone modes, originating from the 
breaking of the translational and the global U(l) phase 
symmetry if one assumes a definite value of the displace- 
ment q, and, e.g., a real order parameter $(x;g). The 
zero modes together with their adjoint mode functions 
are [T3 . (No is the mean atom number in $) 

u q (x;q) — —i K^/n sech [k(x — q)] 
u s (x;q) = $(x;q) 



uf(x;q) 



4V" 

<&(x; q) + ixu q (x; q) 
2 (N + n/n) : 



(10) 



with the associated functions v a (x) = —u^(x) and 
v a d ( x ) ~ u ^*( x )- The quantization box is much larger 
than £, and we neglect small boundary corrections. 

We are now in position to compute the mean atomic 
density and the covariance matrix C in the Bogoliubov 
approximation. The total matter-field operator is split 
into 4f(x) = ®(x) + 8~&(x), with [SV (x) , 5& (y)] = S(x - 
y), and where S$(x) = J2a=6,q PaK d ( x ) ~ iO a u a (x) + 
J2k huk(x) + b\vl(x) with [6 a ,Pp] = i 6^ as in Ref.fli], 
and [6fc, 6+] = 5^. 

The Hamiltonian of the quantized field theory, ex- 
panded to second order in 5^> and , contains 
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quadratic terms P 2 /2m a with effective masses m a , and 
the usual phonon number operator contribution. There- 
fore, in the lowest-energy state (and in thermally excited 
states) of the system, there are no correlations between 
the phonon operators and the zero-mode operators. 
The average density of atoms is given by 

{&(x)*(x)) = \<S>(x)\ 2 + \Mx)\ 2 + Z{x), (11) 



where 

Z(X) = {\ U a(^\ 2 {P") + K(x)\ 2 {9 2 a ) 



(12) 



- u* a (x) < d (x) - 2 3[<(x) uf(x)] (P a § a )} , 

and the dependence on q has not been displayed explic- 
itly. Given the modes (fTU|) , for real <E>, the last term 
in (|12[) vanishes, and the zero-mode contribution to the 
density at the soliton position becomes 



4 (N + n/nf 



{ -^+nn 2 (6 2 )--. 
16 n v ql 4 



(13) 



a sensitivity that scales as n -3 / 4 . In FigQTc), the 
dashed line represents the information that can be ex- 
tracted from the signal-to-noise ratio for a gain func- 
tion g(x,n) = {l — tanh 2 (Ka;)[l + /?(n)]/2} / tanh(ra), 
where (3(n) — tanh(0.014n£ — 0.84) represents a small 
correction to the optimal weighting function for a mean- 
field condensate. 

Why docs the inclusion of noise in the Bogoli- 
ubov description provide an even better resolution than 
the mean field? To understand this, consider as in 
Ref. [ljj the mean value of the detected signal S ~ 
q J dxg(x)d q p(x; 0) (for small g), and its variance AS* 2 = 
J dxdyg(x)g(y)C(x, y\ 0). The amount of information 
that can be extracted from this estimation strategy is 
related to the signal-to-noise ratio S 2 / AS 2 . It can be 
shown from Eq. (fl"4|) . using the completeness relation for 
phonon and Goldstone modes, that AS* 2 = AS^ F + 
AS 2 h — ASq for any gain funtion g(x). The Goldstone 
and phonon contributions, ASq and AS 2 h , are of the 
same order in n, but AS 2 , dominates and we get smaller 
noise 1191. 



IV. CONCLUSIONS 



The operator Pg reflects the fluctuation in the total con- 
densed atom number. Poisson statistics thus implies 
(Pg) ~ N . The soliton position is described by the 

operator 6 q , and behaves like a free particle. We adopt 
here the strategy, introduced in Ref. 20] for a trapped 
BEC, choosing the quantum state for the q mode to be 
the unique Gaussian state in P q , 6 q that minimizes the 
average density n(q) = (^/^ (q)^/ (q)) . This state has a 
finite spread (P 2 ), so it evolves slowly in time, and we 
consider measurements done quickly, before the spread- 
ing of the soliton ultimately invalidates the Bogoliubov 
approach. 

The density correlation function C(x,y;q) is found by 
a straightforward expansion of field operator products to 
second order in <5^, 

C(x,y;q) = $(xWy){J2 

[u k (x) + v k (x)] [ul(y) + v^(y)] 



(P 2 )uf(x)uf*(y) + (9 2 q )u q (x)u q (y) }. (14) 



We finally get the elements of the covariance matrix C by 
integrating C(x,y;q) over finite-size pixels, x £ [xj,Xj + 
Ax] and y £ [x s , x s + Ax]. 

The results of our numerical analysis are presented in 
Figfljc), which shows the Fisher information for both 
the mean-field and the Bogoliubov descriptions of the 
condensate, where the latter is computed within the 
Gaussian approximation. All continuous lines show a 
linear dependence, which implies in our dimensionlcss 
units that the soliton position q can be probed with 



In conclusion, we have used signal processing theory to 
show that the position of a dark soliton can be measured 
with a precision that scales more favorably than the usual 
shot noise. Our study thus illustrates that atomic inter- 
actions can enhance the performance of atom interferom- 
eters. Scaling below the shot-noise limit is also reported 
for entangled or squeezed states, but such states are not 
required in our scheme. Note that a similar potential 
for high-precision sensors may be achieved with bright 
solitons [21( and with vortex lattices, or with a "bubble" 
floating in a BEC and filled with another atomic species 
[H. The n~ 3 / 4 scaling can be understood by simple sta- 
tistical arguments applied to the mean-field condensate. 
Our quantitative analysis by means of the Cramer- Rao 
bound shows, however, that quantum fluctuations due to 
phonon and Goldstone modes enhance the sensitivity be- 
cause these modes provide density correlations that are 
sensitive to the soliton position and make a larger amount 
of information available. 
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